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Abstract 

To understand the mechanisms by which the re-solution of Fe and Cr additions increase the corrosion rate of irradiated 
Zr alloys, the solubility and clustering of Fe and Cr in model binary Zr alloys was investigated using a combination 
of experimental and modelling techniques — atom probe tomography (ART), x-ray diffraction (XRD), thermoelectric 
power (TER) and density functional theory (DFT). Cr occupies both interstitial and substitutional sites in the a-Zx 
lattice; Fe favours interstitial sites, and a low-symmetry site that was not previously modelled is found to be the most 
favourable for Fe. Fattice expansion as a function of Fe and Cr content in the a-Zx matrix deviates from Vegard’s 
law and is strongly anisotropic for Fe additions, expanding the c-axis while contracting the f 2 -axis. Matrix content 
of solutes cannot be reliably estimated from lattice parameter measurements, instead a combination of TER and ART 
was employed. Defect clusters form at higher solution concentrations, which induce a smaller lattice strain compared 
to the dilute defects. In the presence of a Zr vacancy, all two-atom clusters are more soluble than individual point 
defects and as many as four Fe or three Cr atoms could be accommodated in a single Zr vacancy. The Zr vacancy is 
critical for the increased apparent solubility of defect clusters; the implications for irradiation induced microstructure 
changes in Zr alloys are discussed. 


1. Introduction 

Zr alloys are widely used in the nuclear industry as 
nuclear fuel cladding and other structural components. 
Fe and Cr are common alloying elements, added to im¬ 
prove corrosion resistance CHU. These elements are 
known to exhibit near-negligible solid solubility in a- 
Zr, and therefore segregate to form second phase par- 
tides (SPPs) mil. One of the key aspects of the mi- 
crostructural evolution of Zr alloys under irradiation is 
the amorphisation and subsequent dissolution of SRRs, 
which leads to the re-solution of the alloying elements 
in the Zr matrix above their solubility limits 11711191 . In 
turn, this has an impact on the physical and corrosion 
properties of the alloy. In particular, surface oxidation 
and hydrogen pick-up fraction are known to be strongly 
affected by alloy composition and the presence and dis¬ 
tribution of SRRs and experience a marked increase after 
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SRR dissolution |[ni3l [20ll26l . It is important to limit hy¬ 
drogen uptake during reactor operation because hydro¬ 
gen causes dimensional changes to the cladding ED, 
reduces its ductility l24l and reduces integrity perfor¬ 
mance in hypothetical accident scenarios ||27l[28l, and 
potentially in the storage conditions relevant to spent 
nuclear fuel 17711291 . 

Recent advanced transmission electron microscopy 
(TEM) and atom probe tomography (ART) stud¬ 
ies ED have shown that clusters of Fe and Cr form at 
(a) and (c) dislocation loops following the re-solution 
process. This was previously suggested by TEM inves¬ 
tigation ||9 UTT1IT7|[T^I^[^ but not observed directly. 
It has been suggested that irradiation induced defects 
may also act as trapping sites for hydrogen, thereby in¬ 
creasing the terminal solid solubility of hydrogen in a- 
Zr |[34l|35l. 

The solubility of Fe in Zr — and to a lesser extent 
also that of Cr in Zr — has also been investigated us¬ 
ing atomic scale simulations, but so far, the clustering 
behaviour of alloying elements has hardly been con- 
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sidered using such methods. Early work by Perez and 
Weissmann 1361 and of Pasianot et al (371 investigated 
the possible mechanism for Fe accommodation in the 
a-Zx lattice, but their DFT calculations were limited to 
small supercells containing 36 and 48 Zr atoms respec¬ 
tively. In particular, Pasianot et al. Ell observed that 
when Fe substitutes for Zr, it occupies a low symmetry 
configuration that is displaced slightly from the lattice 
site. More recent studies have employed, in one case, 
slightly larger supercells (54 Zr atoms by Fumley et al. 
151 and 48 Zr atoms by Christensen et al. (381133), but 
only the more conventional interstitial sites (tetrahedral 
and octahedral) were considered. Furthermore, recent 
publications l40H46l have shown that even larger super¬ 
cells (~ 300 atoms if no finite size correction term is ap¬ 
plied) are required to avoid computational artefacts that 
may significantly affect the apparent stability of defects 
in a-Zx. It is evident that a state-of-the-art evaluation of 
the extrinsic defects in Zr is needed. 

In previous papers, the authors considered the forma¬ 
tion of SPPs (6l[43 and their interaction with H l^l48l . 
Here, the authors are concerned with the conditions rel¬ 
evant to irradiated Zr alloys, in which the SPPs are par¬ 
tially dissolved. The current work employs a suite of 
experimental and theoretical approaches to investigate 
the solubility of Fe and Cr in pristine and defective Zr 
and the formation of clusters containing Fe, Cr and in¬ 
trinsic defects. First, DFT simulations reveal that Cr 
may occupy both interstitial and substitutional sites in 
the a-Zx lattice, and the results were corroborated by 
spatial distribution maps produced with APT. The sim¬ 
ulations also indicated that the Fe-Zr binary systems ex¬ 
hibits a large deviation from Vegard’s law, thereby indi¬ 
cating that lattice parameter measurements by XRD do 
not provide a suitable estimate of solute concentration 
in the a-Zx matrix. The matrix content of fast quenched 
binary samples was then measured using TEP and APT, 
showing that an increasing amount of Fe and Cr was 
trapped in solution with increasing nominal composi¬ 
tion, despite the formation of SPPs. The lattice ex¬ 
pansion due to alloying additions of the binary samples 
was then measured by XRD, and the predicted devia¬ 
tion from Vegard’s law was observed. Complementary 
DFT simulations also highlight that the preference for 
interstitial over substitutional accommodation is depen¬ 
dent on the atomic strain environment and an argument 
is put forward for clustering of dilute defects as a means 
to reduce overall lattice strain. A first nearest neigh¬ 
bour analysis carried out with APT provides evidence 
that these clustering tendencies occur in oversaturated 
Cr-Zr alloys. Further simulations show that larger de¬ 
fect clusters may be favourably accommodated in the 


a-Zx lattice and that Zr vacancies are crucial for the for¬ 
mation and growth of the clusters. Finally, the work is 
summarised and the implications for irradiated Zr alloys 
discussed. 

2. Methodology 

2.1. Materials 

Binary Zr-alloys were melted in an arc furnace in 
a water cooled copper crucible under an argon atmo¬ 
sphere at Western Zirconium, USA. The Zr starting ma¬ 
terial was in the form of chips while Cr and Fe were 
small beads. All alloying elements were standard ma¬ 
terials used by Western Zirconium for their production 
of zirconium alloys. The 125 g buttons were re-melted 
three times in order to ensure chemical homogeneity. 
Further details of the materials processing can be found 
in (43|50]|. The buttons were analysed at Western Zir¬ 
conium using induced coupled plasma-atomic emission 
spectroscopy and combustion analysis (oxygen and ni¬ 
trogen) to determine the chemical constituent of each 
sample, presented in Table It is acknowledged that 
sample Zr-0.05Cr contains notable amounts of Fe and 
Sn contaminations and the results from this alloy are 
highlighted in subsequent section. All other alloys were 
produced with a high degree of purity. 

Table 1: Chemical composition of the binary alloys in 
wt. ppm. Hf and Nb were consistently less then 23 and 
20 ppm respectively. Si and any other potential impuri¬ 
ties were always below the detection limit. 


Sample name 

Cr 

Fe 

Cu 

N 

0 

Sn 

Zr-O.lFe 

<20 

1049 

10 

NA 

NA 

< 8 

Zr-0.2Fe 

<20 

1927 

11 

NA 

NA 

< 8 

Zr-0.4Fe 

< 20 

4298 

< 10 

44 

810 

< 8 

Zr-0.6Fe 

<20 

6226 

22 

NA 

NA 

< 8 

Zr-0.8Fe 

<20 

8943 

19 

NA 

NA 

< 8 

Zr-0.05Cr 

475 

217 

11 

NA 

NA 

1155 

Zr-0.15Cr 

1608 

37 

< 10 

NA 

NA 

< 8 

Zr-0.30Cr 

2869 

41 

< 10 

43 

849 

< 8 


The as-cast buttons were cross rolled at 540 °C with 
an intermediate recrystallisation anneal at 600 °C to a 
final thickness of 3 mm. Subsequently, 3 x 3 x 40 
mm^ matchstick samples were cut and p heat-treated 
for 10 minutes at 1000 °C in a vertical furnace flushed 
with argon, followed by water quenching, in an attempt 
to maintain most of the Fe and Cr into a-Zx solution. 
Scanning and transmission electron microscopy inves¬ 
tigation showed that complete solid solutions were not 
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obtained even at these very high cooling rates. Instead 
a significant number of small SPPs had formed ll50l . 

2.2. X-ray diffraction 

The matchsticks were cut into 3x3x2 mm^ cubes, 
mounted in a 5 x 5 grid before grinding and polishing to 
produce an XRD-sample with an approximate surface 
dimension of 15 x 15 mm^. The x-ray analysis was car¬ 
ried out on a Philips PW3710 diffractometer using Cu- 
Ka radiation in Bragg-Brentano geometry. Diffraction 
profiles were recorded ranging from the {lOTO} to the 
{3032} reflections with a step size of 26 = 0.02° and a 
recording time of 20 s per step. The profiles were anal¬ 
ysed using Rietveld analysis to determine the a and c 
lattice parameters. 


2.3. Thermoelectric power measurements 

TEP experiments measure the Seebeck coefficient 
(S), which is the electric potential difference that arises 
when two metals in tight contact form a thermocouple 
with two junctions held under a temperature difference. 
The Seebeck effect consists of two parts: a chemical 
gradient found at the junctions between the two met¬ 
als (the Peltier effect) and a thermal gradient within 
the same metal (the Thomson effect). In the present 
case, matchstick samples were clamped between cop¬ 
per blocks, which were maintained at temperature T and 
T -I- AT, respectively. The measured TEP, which is rel¬ 
ative to the reference metal to which it is clamped, can 
then be plotted as a function of concentration of solid 
solution in the matrix ETHH. The sign of the TEP can 
be discussed in terms of incomplete d-bands of electrons 
in transition metals m and Eermi surfaces 1561 . Eor an 
element with unfilled d-orbitals, such as Zr, the addition 
of a solute atom into the matrix can either increase or 
decrease the TEP of the alloy. As all transition met¬ 
als have fewer available d-orbitals than Zr, the energy 
difference decreases, which results in the TEP becom¬ 
ing more negative with increasing solute concentrations 
(571. Previous work has shown that the Seebeck co¬ 
efficient of zirconium alloys is sensitive to solute con¬ 
centration, texture and cold work, but is not affected by 
the presence of small volume fraction of SPPs ll58H^ . 
All samples were prepared using the same procedure to 
minimise variations in texture and micro structure (see 
section [TT]), therefore any change recorded in TEP can 


be attributed to variations in solute concentrations. 

The measurements were conducted at INS A, Lyon, 
Erance using a TechMetal Promotion instrument and a 
Cu reference. The temperature of the clamping copper 
blocks was held at 15 + 0.2 °C and 25 + 0.2 °C. In order 


to stabilise the thermo-electricity, each specimen was 
left for 1 min after mounting before measurement. Each 
measurement had a duration of 20 s in which the initial 
value and the variation from this value were recorded. 
Each surface of the matchstick specimen was measured 
twice to give an average value from eight measurements 
per alloy concentration. 

2.4. Atom probe tomography 

APT was performed on a Cameca LEAP 3000 X Si, 
with a flight path of 90 mm. The experiments were con¬ 
ducted at a base temperature of -213 + 5°C, in laser¬ 
pulsing mode (~ 10 ps, 532 nm, spot size < 10 [im di¬ 
ameter). Throughout the analysis, the DC voltage was 
increased to keep a detection rate of 5 ions per 1000 
pulses. Specimens were prepared by means of a EIB 
lift-out procedure, from a mechanically polished sample 
of the Zr-alloys, using a Zeiss Auriga and an electropol- 
ished Mo grid as support (621. The datasets were recon¬ 
structed using state-of-the-art algorithms (631, resulting 
in the typical tomograms shown in Pigure[2a) for a bi¬ 
nary Cr-Zr specimen containing 0.26 at % Cr. Therein, 
a 2 at% Cr isoconcentration surface highlights the pres¬ 
ence of small regions enriched in Cr, up to 4-5 at%, 
which appear to be aligned, maybe along a twin bound¬ 
ary, similar to what was discussed in (49l . The signal 
to background noise ratio on the major peak of Cr was 
120:1, and on the major peak of Ee it was above 200:1, 
providing a high level of certainty when labelling the 
Cr and Ee atoms. Isotopic ratio of ^^Ee/^'^Ee was found 
to be 14.3, which compares well with the natural abun¬ 
dance (15.81). This provides confidence that the signal 
was mostly generated from Ee ions and not molecular 
impurities such as CO. 

As evidenced in Eigure which contains a two- 
dimensional density map computed from the same to¬ 
mogram, the data exhibits features that can be directly 
related to the crystallography of the specimen (64l : a 
so-called pole can be seen at the bottom left along with 
a series of zone lines exhibiting six-fold symmetry. This 
pole can be attributed to the (0002) atomic planes at the 
specimen surface causing trajectory aberrations. These 
planes are also imaged within the APT data. Imag¬ 
ing atomic planes is common in APT, and allows for 
calibration of the tomogram (65l as well as site 
occupancy analyses Elffiol, which are facilitated by 
data treatment methods such as spatial distribution maps 
(^|7T]|72l. The latter are similar to split radial distribu¬ 
tion functions investigating the local neighbourhood of 
each atom along a specific direction. Another common 
data treatment technique, the calculation and analysis of 
the distance between a given species and its first nearest 
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Figure 1: (a) Three-dimensional reconstruction of the dataset from the Zr-Cr sample. The green surface encompasses 
regions containing above 2 at.% Cr. (b) Top-down projection of the dataset shown in (a), with red representing regions 
of higher density and blue of lower density. 


neighbours CSEl, was also used here to estimate the 
matrix composition following the protocol described in 
ca as reported in B9l . but also to investigate the clus¬ 
tering tendency of Cr in Zr. 

2.5. Computational calculations 

DFT simulations were carried out using the castep 
code |[76]| with the PBE exchange-correlation functional 
ca, ultra-soft pseudo potentials fTSll and a consistent 
plane-wave cut-off energy of 450 eV. 

Supercells containing 150 Zr atoms were modelled 
using a 2 X 2 X 2 k-point sampling grid 11791 . The 
linear elastic theory correction term of Varvenne et al. 
Il45l {anew) was employed to reduce finite size effects. 
The elastic constants of a-Zr, which feed into anew, 
were calculated by performing small lattice perturba¬ 
tions. The resulting stiffness constants are (in units 
of GPa): cii = 141.97, ci 2 = 65.36, cb = 68.02, 
C 33 = 148.71, C 44 = 30.22 and C 66 = 38.30. Since these 
systems are metallic, density mixing and Methfessel- 
Paxton 1801 cold smearing of bands were used (smear¬ 
ing width = 0.1 eV). Tests were carried out to ensure a 
convergence of 1 x 10“^ eV/atom with respect to all pa¬ 
rameters. No symmetry operations were enforced when 
calculating point defects and all simulations were spin 
polarised. 

The energy convergence criterion for self-consistent 
calculations was set to 1 x 10“^ eV. Similarly robust 


criteria were imposed for atomic relaxation within the 
memory conservative BFGS algorithm l8TJ[82l: the en¬ 
ergy difference was less than 1 x 10“^ eV, forces on in¬ 
dividual atoms less than 0.1 eV nm“^ and for constant 
pressure calculations, stress components on the cell of 
less than 1 MPa. 

Defect formation energies were calculated using 
equation [2 

Ef = £DFT _ £DFT + ^ ^ 1 

i 

where and E^^ are the total energies from the de¬ 
fective and perfect DFT simulations, pt is the chemical 
potential of all species i that are added or removed from 
the perfect crystal to form the defect, and £’int is the cor¬ 
rection term for the interaction energy of the defect with 
its periodic images, calculated using anew Ea. The 
chemical potential p is calculated as the DFT energy 
per atom of the metallic elements in their ground state; 
for Fe the ground state is the ferromagnetic BCC phase, 
for Cr it is the anti-ferromagnetic BCC phase. Anew 
calculations employed a radial cutoff of 15 A with 20 
divisions of the Fourier grid, which yielded energy val¬ 
ues converged up to the 4* decimal place. 

The relaxation volume (Q) of a defect is defined as 
the difference in volume between a supercell containing 
the defect and the perfect Zr supercell; see equation 

Q = V{Zr,My) - V{Zr,) (2) 
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When calculating Q, mass action is not taken into ac¬ 
count, in other words, the number and types of atoms 
between the defective and perfect cell do not have to 
be the same. In fact, considering the subscript of equa¬ 
tion for an interstitial defect v = z, whilst for a sub¬ 
stitutional defect X y = z. A related quantity often 
found in the literature is the defect formation volume, 
in which the number of atoms of each species is con¬ 
served between perfect and defective cells. However, 
the defect formation volume is only properly defined for 
intrinsic defects (831, as the reference volume of an iso¬ 
lated extrinsic atom is not a strictly defined quantity. In 
the current work, only relaxation volumes will be con¬ 
sidered. 

Configurational averaging of physical properties re¬ 
lated to the presence of defects (such as Q) was per¬ 
formed following equation!^ 


2 ^ Z/^(0^(0exp(-a) 
2/^(0 exp(-2,) 


(3) 


where 

_ AEfii) 
keT 


(4) 


and X is the physical property of interest, i is a defect 
configuration, n(i) is its multiplicity, AEf(i) is the differ¬ 
ence in formation energy with respect to the most stable 
defect and all other symbols retain their conventional 
meaning. 


3. Results and Discussion 

3.1. Fe-Zr and Cr-Zr binary series 

TEP measurements of the samples evaluate the See- 
beck coefficient of the a-Zr matrix, bypassing the resis¬ 
tance contributions of the SPPs There- 

fore, any change in Seebeck coefficient observed across 
the binary alloys, compared to the unalloyed reference 
sample, can be related to the amount of Fe or Cr in 
solution. TEP results are plotted against composition 
in Figure It is observed that the Seebeck coefficient 
decreases with increasing content of Fe or Cr, indicat¬ 
ing that, despite the formation of SPPs, an increasing 
amount of alloying element was trapped in solution with 
increasing nominal composition of the samples. 

To obtain quantitative concentrations for the a-Zr ma¬ 
trix, it is necessary to establish a datum by calibrating at 
least one TEP point that has a known alloy concentra¬ 
tion in solution. Since all of the samples exhibited some 
segregation of Fe/Cr to SPPs, APT was employed to cal¬ 
culate the matrix content of Fe and Cr in a-Zv by sam¬ 
pling a volume containing no SPPs. The matrix content. 
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Figure 2: TEP measurements against nominal compo¬ 
sition of the binary alloys. Right-hand side axes are the 
calibrated concentrations of Cr and Fe dissolved in solu¬ 
tion. Error bars express the standard deviation of eight 
samples per datum. The hollow symbol represent the 
low purity alloy. 


therefore, includes all the Fe and Cr atoms that are not 
within SPPs, irrespective of the atoms being uniformly 
dispersed in solution, segregated at the grain boundaries 
or clustered in atmospheres of higher Cr or Fe density. 
These results were published previously (491 , and the 
relevant findings are summarised thus. The samples 
containing 1.30 at.% Fe and 0.26 at.% Cr were chosen 
for APT analysis. Fe and Cr are either found within 
small particles or atmospheres along grain boundaries, 
or are randomly distributed within the matrix. The to¬ 
tal composition calculated using APT is in excellent 
agreement with the nominal composition of the alloys 
(1.34 + 0.026 at.% Fe and 0.25 + 0.01 at.% Cr respec¬ 
tively). Selecting volumes that contained no SPPs or 
grain boundaries, the concentration of alloying element 
in solution was calculated to be 0.42 + 0.015 at.% Fe 
and 0.21 + 0.01 at.% Cr. 

The change in Seebeck coefficient measured by TEP 
is reasonably linear with increased alloy concentration 
(see Figure [^. Therefore, assigning the concentration 
values obtained by APT to the points at 1.30 at.% Fe 
and 0.26 at.% Cr, and extrapolating linearly so that the 
reference sample has an alloy concentration of zero, the 
solution concentration of the remaining samples was 
obtained. The resulting scales are plotted on the right- 
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hand y-axis in Figure The deviation from the linear 
regression then represents a new estimate of the uncer¬ 
tainty, which is added to the uncertainty arising from 
ATP measurements and that of TEP measurements, to 
form the total error about the solution concentration of 
each sample (used in later graphs). In the following sec¬ 
tions, when referring to alloy composition, we refer to 
these calculated concentrations of Cr or Fe in solution 
rather than the nominal composition of the buttons. 

3.2. Fe and Cr solubility in a-Zr 

To calculate the solution energy of the alloying el¬ 
ement within the non-interacting regime, DFT simula¬ 
tions were performed with single Fe and Cr defects in 
the supercells described previously, equivalent to an al¬ 
loy concentration of 0.67 at % Fe/Cr. A number of re¬ 
cent studies Il40ti46l have highlighted that finite cell size 
effects may significantly affect the predicted stability of 
Zr self interstitials atoms (SIAs); simulations performed 
using a small size of supercell, or without the use of 
finite-size correction methods, yield spurious results. In 
particular, the work by Varvenne et al. EH showed that 
with an energy correction term calculated from linear 
elastic theory, as used in the current work (£’/„?), su¬ 
percells containing 200 Zr atoms accurately described 
Zr SIAs, and supercells containing only 96 Zr atoms 
yielded differences of only 40-150 meV. All the defects 
considered in the current work cause significantly less 
lattice strain than Zr SIAs. This provides confidence 
that the supercell employed in this study, which con¬ 
tained 150 Zr atoms, is sufficiently large to avoid spuri¬ 
ous finite size effects. This is corroborated by the fact 
that was in the range of 3-60 meV for all point de¬ 
fects. As a further confirmation, all point defect simu¬ 
lations were repeated under constant volume conditions 
(s = 0) and constant pressure conditions (cr = 0), and 
the difference in energy between the two methods was 
consistently less then 0.6 %. 

Many interstitial positions, as well as the substi¬ 
tutional and Zr-Fe and Zr-Cr dumbbell configura¬ 
tions were considered; The resulting formation ener¬ 
gies relaxation volumes (Q^"^) and anisotropic 

strains on the supercell reported in Ta¬ 

ble 1^ for all stable defectsjj Some interstitial positions 
were found to be unstable, that is, the defects moved 
to another site upon relaxation; these include the tetra¬ 
hedral positions (which appears as stable when simu¬ 
lated in small supercells), the hexahedral position (also 


and were within 0.07% and 3.7% of and 
respectively. 


termed basal tetrahedral by some authors) and all the 
dumbbell configurations. The substitutional defect con¬ 
sistently relaxed to the off-site substitutional position 
discussed in unless symmetry constraints were im¬ 
posed. This suggests that the high spin high-symmetry 
substitutional site observed by Christensen et al. |[38]| 
may be due to insufficiently high degree of convergence 
during geometry relaxation. 

Similarly, the tetrahedral site relaxed into the newly 
observed crowdion configuration if a suitably large sim¬ 
ulation cell is employed. The current work identifies an¬ 
other interstitial site that has not previously been simu¬ 
lated: the off-site octahedral. This is significantly more 
stable than any other interstitial site for Fe, but is un¬ 
stable for the accommodation of Cr, which is consistent 
with the larger atomic radius of Cr. Mossbauer studies 
of Fe in a-Zr (841 suggested that ~ 30% of the total Fe 
in solution is located in off-centre interstitial sites of this 
type. 

Regarding the accommodation of Fe, all interstitial 
sites provide more favourable solution energies than the 
substitutional site. This is in agreement with experimen¬ 
tal diffusivity measurements (85^871 . and all previous 
DFT calculations |[6| |36l [371 with the sole exception of 
(38l . It is unclear from the literature whether Cr atoms 
exhibit a similar preference for interstitial accommoda¬ 
tion. Experimental diffusivity measurements indicated 
that Cr diffuses 2-4 orders of magnitude slower than 
Fe in a-Zr. However, Fe is reported to diffuse 6-9 or¬ 
ders of magnitude faster than for Zr self-diffusion and 
9-16 orders of magnitude faster than substitutional so¬ 
lutes |[8M89l . This suggests that the transport of Cr in 
the a-Zr lattice may be mediated by an interstitial so¬ 
lute. The current work is in excellent agreement with 
a previous DFT publication HI, which highlights that 
whilst the preferred site for Cr solution is substitutional, 
the difference in energy between that and the intersti¬ 
tial octahedral site is very small, and therefore Cr is 
expected to exhibit both substitutional and interstitial 
behaviour. Despite this, subsequent modelling studies 
have not considered the possibility of interstitial accom¬ 
modation for Cr (381. 

Experimental evidence of the dual nature of Cr ac¬ 
commodation in Zr is provided by our APT work. The 
dataset of the sample is shown in Figure ^sl). A 
10 X 10 nm^ subset of the data centred on the (0002) 
pole indicated in Figure [Jb), and going down the whole 
length of the dataset, was exportecQ Advanced species- 
specific spatial distribution maps were applied to this 


^Due to the small size of the region, a large portion of those atoms 
reside at the boundary and have a limited number of neighbours. 
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Table 2: Defect formation energy {Ey^) and volumetric properties for all defects that may accommodate Fe or Cr 
in bulk Qf-Zr. and are the strain in the a and c direction respectively. Full relaxation volume tensors are 


presented in Appendix A 




(eV) 

(A^) 

(%) 

(%) 

Fe 

off-site substitution 

1.388 

-10.40 

-0.15 

-0.15 


substitution 

1.709 

-16.25 

-0.02 

-0.40 


octahedral 

1.079 

13.70 

-0.02 

0.27 


off-site octahedral 

0.941 

13.53 

-0.07 

0.25 


trigonal 

1.212 

13.42 

0.18 

-0.26 


crowdion 

1.172 

13.52 

0.10 

0.09 

Cr 

off-site substitution 

1.732 

-11.31 

-0.03 

-0.21 


substitution 

1.892 

-12.83 

-0.11 

-0.13 


octahedral 

1.882 

15.20 

-0.05 

0.30 


trigonal 

1.968 

13.65 

0.16 

-0.17 


crowdion 

2.061 

15.41 

0.15 

0.08 


subset and the resulting data are plotted in Figurej^ The 
Zr-Zr distribution exhibits broad peaks corresponding 
to the (0002) atomic planes. The slight and progres¬ 
sive shift away from the expected location of the peaks 
can be attributed to distortions in the tomogram, as dis¬ 
cussed in ll90l . The Zr-Cr distribution, which measures 
the average distribution of Cr atoms relative to Zr atoms 
along this crystallographic direction, exhibits peaks that 
are in-between the main peaks of the Zr-Zr distribution. 
This is consistent with a significant fraction of the Cr 
atoms being located at interstitial sites. This is the first 
evidence of interstitials provided by APT. A similar pro¬ 
cedure was attempted on the Fe-containing sample, but 
the large volume fraction of SPPs in the sample m 
made it extremely challenging to achieve a sufficient 
signal-to-noise ratio to generate definitive results. 

3.3. Lattice expansion of binary alloys 

Accommodation of Fe and Cr point defects in the a- 
Zr matrix is predicted to cause noticeable lattice strain 
(see Table [^. Therefore, the lattice parameters of the 
a-Zr matrix were measured for all binary samples using 
XRD. The change in lattice parameter a as a function of 
alloy concentration is shown in Figure together with 
a projection of Vegard’s law calculated from the lattice 
parameters of the pure elements (dashed line) and the 
DFT predictions of lattice expansion due to the pres¬ 
ence of defects (dotted lines). Note that the change in 
composition and the change in lattice parameter are very 
small, near the accuracy limit of the XRD equipment, as 
highlighted by the large error bars. 


Therefore, the results taken do not bare a quantitative weight, but can 
provide qualitative information. 



-1.0 -0.5 0.0 0.5 1.0 

z offset (nm) 


Figure 3: In-depth spatial distribution map showing the 
Zr-Zr distribution (in purple) and the Cr-Zr distribution 
(in blue). Dashed lines indicate J(ooo 2 ) spacings from 
XRD data. 


The DFT predictions rely on configurational averages 
of the dilute point defects at temperatures of 25 K, 300 K 
and 600 K. This averaging technique does not include 
other temperature effects such as thermal expansion or 
phonon scattering. In other words, the model represents 
a solution that has been homogenised at those tempera¬ 
tures and subsequently quenched. 

Figure shows that, for Cr-Zr solution, both Veg¬ 
ard’s law and the DFT predictions are in good agree- 
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Cr at. % Fe at. % 

Figure 4: Measured lattice parameter a as a function of alloying element in solution (points) and lattice parameters a 
and c from theory (lines). Dashed line represents Vegard’s law; Dotted lines represent the current prediction of lattice 
expansion using configurational averaging of DFT simulations of dilute defects. Hollow symbol represent low purity 
sample. Vertical error bars represent the uncertainty due to machine error, background noise and Rietveld analysis; 
horizontal error bars represent the compound uncertainty of TEP measurements, APT concentration measurements 
and standard deviation from linear regression of Figure]^ 
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ment with experimental observations, if the low purity 
sample (hollow circle) is disregarded. However, for the 
Fe-Zr solution DFT predictions differ greatly to Veg- 
ard’s law. Further, for the a lattice parameter, DFT pre¬ 
dictions are in better agreement with experimental data 
at low concentrations (near the solid solubility of Fe in 
Zr) but the agreement is somewhat lost at higher con¬ 
centrations. For the c lattice parameter, DFT results are 
in stark contrast to Vegard’s law in that an expansion 
is expected instead of a contraction. XRD data for the 
c lattice parameter were inconclusive as the low multi¬ 
plicity of the c direction caused too much scatter in the 
data. The predicted deviation from Vegard’s law sug¬ 
gests that lattice parameter measurements are not a suit¬ 
able means to estimate the solute concentration of alloy¬ 
ing elements. 

When performing the configurational average, it is 
implicit that the defects are not interacting; therefore, 
strictly, the average is only valid at the dilute limit. 
Since the binary solid solutions investigated here are 
above their respective solubility limits, it is pertinent to 
assume that the alloying atoms are interacting with each 
other. More specifically, the compressive strain field of 
an interstitial defect is likely to increase the formation 
energy of another interstitial defect in its vicinity, whilst 
reducing that of a substitutional defect (which, due to its 
negative relaxation volume, has a tensile strain field as¬ 
sociated with it). This hypothesis was corroborated by 
repeating the defect simulations in pre-strained super¬ 
cells, see Figure Under a compressive strain, the sta¬ 
bility of substitutional defects (filled squares) increases 
while that of interstitial defects (hollow symbols) de¬ 
creases; and the opposite is true under a tensile strain. 
This helps explain the lack of preference between in¬ 
terstitial and substitutional accommodation that is ob¬ 
served in the APT spatial distribution map of Cr (Fig¬ 
ure [^. 

As well as affecting the relative solution energies, the 
strain fields of the point defects may provide a driv¬ 
ing force for diffusion: defects with opposing strain 
fields may attract each-other at distances of up to a few 
angstroms, whilst defects with same-sign strain fields 
will repel one another. When combined with the ex¬ 
treme mobilities of Fe and Cr (SIIEIIISD, this may lead 
to the formation of defect clusters with a reduced over¬ 
all lattice strain, hence a reduced lattice expansion and 
more favourable solution energy. 

3.4. Cluster formation 

To investigate the formation of Fe and Cr clusters, 
simulations containing two extrinsic species were first 
considered. The starting positions for the clusters were 


Applied pressure (GPa) 
2.5 0 


-2.5 



■ Crzr 

A Crt, 

O CTcrow 

■ Fezr 

A Fetri 

FGoffoct 
O FSoct 
O FGcrow 


Hydrostatic strain (%) 


Figure 5: Energy of solution of Fe (beige) and Cr (blue) 
accommodated as interstitial species (hollow symbols) 
and substitution species (filled squares) as a function of 
hydrostatic strain. The simulation cells were strained 
prior to adding the defect, by applying an external hy¬ 
drostatic pressure, displayed in the secondary v-axis 
above (positive = compressive). 


defined by combining a substitutional defect (Mzr) with 
and an octahedral or off-site octahedral interstitial de¬ 
fect (M/(oct)), since these are the most stable defects 
with opposing strain fields for Cr and Fe respectively 
(from Section [33] ). All such configurations that could 
fit in a 5 X 5 X 3 supercell of a-Zr (150 Zr atoms) 
were investigated, leading to defect-defect separations 
that range from 2.30 A for the first nearest neighbour 
(Inn) to 6.97 A for the 7nn configuration. When consid¬ 
ering mixed Fe-Cr clusters, the Cr atoms were placed 
in the substitutional sites and the Fe atoms in the in¬ 
terstitial sites, {Crzr : Fe/(off-oct)}, owing to the smaller 
atomic radius of Fe and its preference for interstitial 
sites (see section |3.3| ). The simulations of the clus¬ 
ters were relaxed to a high level of force convergence 
(0.05 eV A“^) and the atomic positions were perturbed 
by small amounts in random directions. Furthermore, to 
provide greater degrees of freedom to the simulations, 
these were repeated under cr = 0 conditions as well as 
£ = 0 conditions. This combination ensures that the 
BFGS minimiser ISTl 1821 is unlikely to trap atoms in 
local energy minima. In other words, the starting posi¬ 
tions are just that, and the extrinsic atoms were expected 
to explore the energy surface until lowest energy config- 
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Figure 6: Brown spheres represent Fe atoms, dark blue spheres represent Cr atoms, turquoise spheres represent Zr 
atoms, translucent spheres represent the initial position of selected atoms. 


urations were found. Figureshows the initial and final 
configurations of some clusters that have moved from 
their original lattice sites. In all cases, the resulting de¬ 
fect is an elongated or extended defect, often involving 
one or more Zr SIAs. Most notably the Inn clusters 
moved into a split substitutional (or dumbbell) around 
a Zr lattice site. An in-depth analysis of the dumbbell 
configurations, including a comparison with the intrin¬ 
sic Zr dumbbells from previous work ||46l [92]|, is pre¬ 
sented in [Appendix B 

A summary of the formation energies, binding ener¬ 
gies and relaxation volumes of all the 2-atom defects — 
in their final relaxed positions — are presented in Ta- 
ble[^ In all cases the dumbbell defect is consistently the 
most stable configuration, independent of the species in¬ 
volved, and the relative preference for the dumbbell is 
as high as 0.5-0.6 eV compared to the next most stable 
configuration. Importantly, all configurations up to the 
^th ^ single lattice jump away from the dumbbell 

configuration. Whilst this provides incomplete informa¬ 
tion about kinetics of cluster formation, it does imply 
that multiple paths exist for migrating extrinsic species 
to reach (and be trapped in) the dumbbell configuration. 

With regard to the lattice expansion, all two-atom 
clusters exhibit relaxation volumes that are significantly 


smaller than those of the single-atom dilute defects. 
Furthermore, the combined defect volume of dilute Fezr 
and Fe/ is 3.13 A^, while that of the bound dumbbell is 
only -2.49 A^. Similarly for Cr, the combined volume 
of dilute defects is 3.81 A^ compared with -1.95 A^ for 
the dumbbell. Finally, in the mixed case (in which Cr 
takes the substitutional site and Fe takes the off-oct in¬ 
terstitial site), the combined volume of dilute defects is 
again greater than that of the dumbbell (2.50 A^ against 
-0.42 A^). This suggests that part of the binding energy 
comes from a reduction of lattice strain. 

Notably, many of the Fe-Fe defect pairs exhibit nega¬ 
tive relaxation volumes, resulting in a tensile strain field, 
despite the addition of one extra atom in the supercell. 
In particular, the most favourable configurations (Inn, 
3nn & 4nn) exhibit tensile strain fields arising from re¬ 
laxation volumes of -2.49 A^, -2.39 A^ and -3.59 A^. 

To investigate if more than two Fe or Cr atoms could 
be accommodated in or around a single Zr vacancy, a 
third and then a fourth interstitial atom were added to 
the relaxed dumbbell configurations, as these are the 
most stable 2-atom clusters (see Figure |7]). The result¬ 
ing solution energies (normalised per extrinsic atom, as 
defined in ll48]| ) and defect volumes are presented in Fig¬ 
ure Clusters containing 3 and 4 Fe atoms exhibit 
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similar solution energies to the 2-atom clusters. Fur¬ 
thermore, from Figure (top) it is clear that even the 
least stable 2-atom clusters exhibit solution energies that 
are similar to the most stable dilute defect. The relax¬ 
ation volumes of the clusters containing a Zr vacancy 
shows a fairly linear behaviour, see Figure [^(bottom). 
As a result, clusters containing 3 Fe or Cr atoms cause a 
smaller expansion of the lattice than a dilute interstitial 
Fe or Cr defect. 


Table 3: Normalised defect formation energy 

(£’//atom), binding energies (£’/,), spin and relaxation 
volumes (Q), for all defect clusters investigated. The 
spin is the cumulative spin on the extrinsic elements cal¬ 
culated using Mulliken analysis 19^ . c{hkjl\ stands for 
a crowdion defect along the [hkjl] direction. Nomen¬ 
clature of dumbbells is defined in ll46l and in [Appendix 


E 



(a) {3M}zr (b) {4M}zr 



cluster type 

spin (h) 

Ef/atom (eV) 

Eb (eV) 

a(Ab 

Fe-Fe 

Inn 

dumbbell (PS’) 

0.00 

0.32 

1.69 

-2.49 

2nn 

{Fezr : Fe,} 

1.50 

0.79 

0.75 

1.00 

3nn 

c[10l0] 

0.00 

0.62 

1.08 

-2.39 

4nn 

c[1120] 

0.00 

0.57 

1.18 

-3.59 

5nn 

c[2021] 

0.70 

0.72 

0.88 

-1.95 

6nn 

{Fezr : Fe,} 

1.64 

0.93 

0.47 

2.32 

7nn 

{Fezr : Fe,} 

1.65 

0.96 

0.42 

2.00 

Cr-Cr 

Inn 

dumbbell (P2S) 

0.00 

0.95 

1.95 

1.92 

2nn 

{Crzr : Cr,} 

2.30 

1.26 

1.34 

7.13 

3nn 

c[1010] 

0.00 

1.41 

1.03 

1.81 

4nn 

c[1120] 

2.26 

1.31 

1.22 

5.52 

5nn 

c[2021] 

0.00 

1.54 

0.76 

2.07 

6nn 

{Crzr : Cr,} 

2.29 

1.37 

1.11 

2.96 

7nn 

{Crzr : Cr,} 

2.29 

1.36 

1.11 

7.91 

Fe-Cr 

Inn 

dumbbell (PS’) 

0.00 

0.64 

1.77 

-0.42 

2nn 

{Crzr : Cr,} 

0.00 

1.17 

0.70 

0.66 

3nn 

(Crzr : Cr,} 

2.32 

0.89 

1.26 

6.89 

4nn 

c[1120] 

2.14 

0.87 

1.31 

-0.20 

5nn 

c[2021] 

2.19 

0.76 

1.52 

5.40 

6nn 

(Crzr : Cr,} 

0.00 

1.36 

0.32 

2.96 

7nn 

(Crzr : Cr,} 

2.27 

0.90 

1.25 

6.34 


Figure 7: Brown spheres represent Fe or Cr atoms, 
turquoise spheres represent Zr atoms and the translucent 
cube represents the Zr vacancy. 

Clustering behaviour of Cr was also investigated us¬ 
ing APT. Two distinct analyses were performed to in¬ 
vestigate the distribution of distance between each Cr 
atom and their first nearest neighbour Cr: first consider¬ 
ing the whole dataset, and second considering a subset 
of the data that excludes the regions delineated by an 
isoconcentration surface similar to the one displayed in 
Figure with a threshold of 1 at% Cr. The latter allows 
for an analysis of the matrix. The two graphs, shown 
in Figure exhibit a different behaviour: for the com¬ 
plete dataset, there is a clear tendency for neighbours of 
shorter distances compared to random, while Cr in the 
matrix are close to a random distribution, which was ex¬ 
pected based on visual inspection. The closer distances 
in the subset containing the clusters, is in agreement 
with the DFT predictions that closely bound clusters 
(those in Inn configuration) are more stable than clus¬ 
ters in which the Cr atoms are further apart and com¬ 
pared to dilute defects. 

Interestingly, and as shown in Figure[^ the Cr-rich re¬ 
gions are seen to align along specific directions, which 
could be twin or low-angle boundaries (a high-angle 
grain boundary would have been identifiable in the 
dataset via a change in the pole structure |[64l ). or dis¬ 
location loops, which are likely to be defective areas 
within the material. This could imply that for Cr to clus¬ 
ter there needs to be some defects in their vicinity. 
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Figure 8: Solution energy (top) and relaxation volume 
(bottom) as a function of cluster size. Dilute interstitial 
defects are also included for comparison. Joined points 
show the lowest energy clusters. Blue crosses represent 
Cr defects, yellow dots represent Fe defects. 


3.5. Implication for re-solution ofFe and Cr in irradi¬ 
ated Zr 

A marked increase in corrosion rates of Zr cladding is 
observed at high fuel burn-up (Him. Due to the delay 
of the irradiation effect, the increase in corrosion rate is 
attributed to the amorphisation and dissolution of SPPs 
that also occurs at high burnup mM- Here we consider 
the fate of Fe and Cr, that are released into the a-Zr 
matrix following SPP dissolution. 

DFT and APT results provide strong indications that 
in the presence of a Zr vacancy, the formation of clusters 
significantly lowers the solution energy of Fe and Cr in 
a-Zr, thereby increasing their apparent solubility in a- 
Zr. The extremely favourable binding energies observed 
for all the two-atom clusters suggest that if a cluster is 



0 1 2 3 4 5 

Nearest neighbour 


Figure 9: Cr first nearest neighbour distribution (solid 
lines) compared to randomised datasets (dotted lines) 
for the whole dataset (blue) and for the matrix only (pur¬ 
ple). Inset is the difference between the experimental 
and random that highlights the clustering tendency of 
Cr in the full dataset. 


formed, there will be a barrier against separation. If the 
clusters do not separate, the only other fate is growth. 
It is well known that Fe and Cr are ultra-fast diffusers 
in a-Zr (ST] [89l |9T1, it is, therefore, likely that a third 
atom will fall within the proximity of the cluster and 
also get trapped. Binding energies were found to remain 
strongly positive with the addition of third or fourth Fe 
atoms. Whilst the simulation of larger clusters is unfea¬ 
sible due to computational limitations, it is reasonable 
to expect that larger clusters may be favourable if more 
that one Zr vacancy was present (e.g. considering the 
interaction of Fe or Cr with the vacancy clusters inves¬ 
tigated by 1941 , or with gliding dislocations MM)- 
Since the mobility of Fe and Cr is not likely to be 
a limiting step in cluster formation, the most probable 
limiting factor is the presence of Zr vacancies. Fe in 
solution in Zr naturally takes an interstitial position (see 
sectionjr^, whilst all the clusters considered in the cur¬ 
rent work contain one Zr vacancy. Therefore, for the 
formation of an Fe/Cr cluster, an additional Zr vacancy 
must be present as well as at least two interstitial atoms. 
Further confirmation of this comes from APT, which 
suggests the need for defects nearby, for the formation 
of Cr clusters. 
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Under equilibrium conditions the concentration of 
Zr vacancies is very low (our DFT calculated vacancy 
formation energy of 1.9 eV yields a concentration of 
3.7 X 10“^ at the yS ^ Of transition temperature), and 
therefore only a limited number of defect clusters may 
form. On the other hand, upon irradiation the con¬ 
centration of vacancies and dislocation loops increases 
by orders of magnitude ifTTIl . Neutron irradiation also 
causes SPPs to become amorphous and to leach Fe and 
Cr back into the Of-Zr matrix (TIBinHBlIISl. In view 
of the current results, this re-solution can be explained 
by the presence of additional Zr vacancies introduced 
through irradiation, which promotes the formation of 
defect clusters. From this point of view, cluster forma¬ 
tion may be considered as a competing mechanism to 
SPP formation. The classical (macroscopic) notion of a 
single a-Zr phase with a uniform random distribution of 
alloying additions is not representative of the irradiated 
Zr-Fe-Cr system at an atomic scale. A better descrip¬ 
tion is one consisting of nano-sized clusters containing 
most of the dissolved Fe and Cr, separated by volumes 
of near pure a-Zr. 

4. Summary 

Binary Fe-Zr and Cr-Zr alloys were investigated us¬ 
ing a combination of experimental and modelling tech¬ 
niques. Buttons of binary alloys were manufactured 
with varying concentrations of alloying elements and 
then quenched at fast cooling rates in an attempt to hin¬ 
der nucleation of SPPs. TEM investigation shows that 
SPPs had formed in all samples. Nonetheless, TEP mea¬ 
surements revealed that an increasing amount of the al¬ 
loying elements were trapped into solution with increas¬ 
ing nominal composition of the samples. APT was em¬ 
ployed to measure the solution concentration of one Fe- 
Zr sample and one Cr-Zr sample. A linear regression 
following the trend of the TEP results was used to esti¬ 
mate the solution concentration in the other binary sam¬ 
ples. 

Results from APT demonstrated that Cr can occupy 
interstitial sites, and that Cr tends to cluster along lin¬ 
ear features, which are likely to be related to defective 
regions of the sample. In turn, this provides confidence 
in the DFT results, which predicted very close energies 
for the accommodation of Cr as an isolated interstitial or 
substitutional species, as well as that Cr tends to cluster 
near vacancies. On the other hand, Fe atoms in dilute 
conditions exhibit a distinct preference for only intersti¬ 
tial accommodation. Furthermore, two low-symmetry 
interstices (the crowdion and the off-site octahedral) 


were successfully modelled for the first time. The off¬ 
site octahedral was found to be the most favourable site 
for the accommodation of Fe, whilst the tetrahedral site, 
which has previously been modelled, was found to be 
unstable. 

XRD results show that limited expansion or contrac¬ 
tion of the a lattice parameter is observed with increas¬ 
ing Fe or Cr content. Vegard’s law provides a reason¬ 
able approximation to the change in lattice parameter a 
for Cr. However, for Fe-Zr solutions, a deviation from 
Vegard’s law is observed in agreement with DFT sim¬ 
ulations. The same calculations also indicate a largely 
different behaviour for the c lattice parameter, which is 
predicted to expand whilst Vegard’s law predicts a con¬ 
traction. DFT predictions based on single point defects 
fall short at higher Fe concentrations, where the dilute 
non-interacting assumption is no longer valid. At those 
higher concentrations, the formation of small clusters 
composed of defects with opposing strain fields is pre¬ 
dicted. 

DFT calculations show that 2-atom clusters exhibit 
much smaller defect volumes compared to their dilute 
defect components. A variety of different clusters were 
identified which were strongly bound and more stable 
than their dilute counterparts, thereby reducing the so¬ 
lution energy of Fe and Cr in a-Zr. In particular, dumb¬ 
bell configurations (i.e. two extrinsic elements around a 
Zr vacancy) were most favourable. The nearest neigh¬ 
bour distribution from atom probe microscopy confirms 
a prevalence of dumbbell configurations for Cr clusters. 
The formation of clusters is predicted to be limited by 
the presence of intrinsic Zr vacancies. This helps ex¬ 
plain the re-solution of Fe and Cr in a-Zr, which is ob¬ 
served upon neutron irradiation, as a direct consequence 
of a higher intrinsic defect population. The current work 
predicts that up to four Fe atoms or three Cr atoms can 
be accommodated around a single Zr vacancy. Larger 
clusters involving more than one Zr vacancy — which 
are currently computationally impractical to model — 
are likely to form in the vicinity of strain gradients 
(caused by grain boundaries, dislocations, voids, etc.) 
and could further reduce the solution energies of Fe and 
Cr. 
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Appendix A. Relaxation tensors of point defects 

The relaxation volume tensors of the dilute point de¬ 
fects from DFT simulations were calculated using the 
aneto program and are presented in Table A.4 The val¬ 
ues are taken from £ = 0 simulations. 


The mixed element dumbbell exhibits similar proper¬ 
ties to the Fe-Fe dumbbell: it adopts the PS’ orientation 
with a tilt angle of ^ = 65.9°, but in this case the dumb¬ 
bell is in-plane but slightly off-centre (shifted 0.321 A 
in the c-direction). 

The difference between the dumbbells investigated 
here and SI A dumbbells, are thought to be related to 
the different ionic sizes and electronegativity of the ex¬ 
trinsic elements. Therefore, whilst the differences may 
appear to be subtle, they are important for the devel¬ 
opment of accurate larger scale models, such as classi¬ 
cal potentials for molecular dynamics, or kinetic Monte 
Carlo simulations. 


References 


Appendix B. Dumbbell configurations 

Here we employ the nomenclature introduced by 
Verite et al 1461 , to refer to the orientation of the dumb¬ 
bell with respect to the Zr lattice: vertical (S), basal 
(BS), rotated on prismatic plane type I (PS or PS’ de¬ 
pending on the rotation angle) and rotated on prismatic 
plane type II (P2S). 

The Fe-Fe dumbbell is rotated in the prismatic plane 
type I (1120) by a tilt angle of 54.5° from the c-axis, 
in accordance with the PS’ defect observed in Zr SIAs. 
The extrinsic dumbbell exhibited lower symmetry than 
the intrinsic SIAs: the Fe dumbbell is ~ 0.25A out of 
plane with respect to the (1120) plane, and also off- 
centre with respect to the host Zr vacancy so that the 
one Fe atom is closer to the vacancy than the other (the 
centre of the dumbbell is displaced by ~ 0.005A with 
respect to the host Zr vacancy). These simulations were 
repeated up to five times with different starting config¬ 
urations and tight force convergence criteria to ensure 
that such subtle measurements were not artefacts caused 
by computational parameters. One other stable position 
was observed for the Fe dumbbell when rotated so that it 
would align vertically. Again the dumbbell is preferen¬ 
tially located ~ 0.40A out of plane. This configuration 
is 0.861 eV less stable than the one described above, and 
was found to lie on a shallow minima near the transition 
state for rotation of the dumbbell. Unlike Zr SIAs, no 
other stable orientations were observed. 

The Cr-Cr dumbbell exhibits a configuration similar 
to the P2S defect (as opposed to the PS’ of the Fe-Fe 
dumbbell), with a large tilt angle of ^ = 76.4°. However, 
the defect is not exactly on the type II prismatic plane 
(lOlO), instead it exhibits a small deviation angle ofif/ = 
8.2° from the (lOlO) plane. Furthermore, the centre of 
the dumbbell is 0.172 A below the Zr lattice. 
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Table A.4: Components of the relaxation volume tensor in units of A^. 




Qii 

^22 



fll3 

^23 

Fe 

off-site substitution 

-2.15 

-4.12 

-4.53 

1.71 

0.18 

-0.32 


substitution 

-3.94 

-3.94 

-8.36 

0.00 

0.00 

0.00 


octahedral 

4.69 

2.87 

6.42 

1.58 

-0.08 

0.14 


off-site octahedral 

4.03 

3.84 

5.95 

2.75 

-0.51 

-0.79 


trigonal 

5.76 

6.67 

-0.18 

3.20 

0.00 

0.00 


crowdion 

5.14 

4.52 

3.86 

1.69 

1.71 

2.92 

Cr 

off-site substitution 

-1.90 

-3.85 

-5.57 

-0.16 

0.00 
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